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The effect of the complex phase of the fermion determinant is a key question related to the 
sign problem in finite-density QCD. Recently it has been shown that ignoring the complex phase 
- the phase quenching - does not change physics in a certain region of the phase diagram when a 
number of colors Nc is large. In this paper we study the effect of the phase quenching within the 
frameworks of effective models and holographic models. We show, in a unified manner, that the 
phase quenching gives exact results for any fermionic observables (e.g., chiral condensate) in the 
mean-field approximation and for gauge- invariant gluonic observables (e.g., Polyakov loop) to one- 
meson-loop corrections beyond mean field. We also discuss implications for the lattice simulations 
and confirm good quantitative agreement between our prediction and existing lattice QCD results. 
Therefore the phase quenching provides rather accurate answer already at A'^c = 3 with small 1/Nc 
corrections which can be taken into account by the phase reweighting. 
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I. INTRODUCTION AND SUMMARY 

Phases of matter under extreme conditions, such as 
the hottest matter in the early Universe and relativistic 
heavy ion collisions, and the most dense matter inside the 
core of neutron stars, are described by Quantum Chro- 
modynamics (QCD) at finite temperature and/or finite 
density. Due to the strong-coupling nature of QCD, lat- 
tice simulations based on importance sampling have been 
the main first principle method to reveal the properties of 
such systems. A number of important properties of hot 
QCD matter, such as the equation of state pj] and a rapid 
crossover from hadronic matter to quark matter 0] , have 
been unraveled near zero chemical potential. However, 
studies of dense QCD matter are difficult because of the 
notorious sign problem. 

Recently it has been realized that QCD at finite baryon 
chemical potential /is (QCDs) is equivalent to QCD at 
finite isospin chemical potential m (QCD/) Q with a large 
number of colors Nc in a certain region of the phase dia- 
gram presumably relevant to the heavy ion collision ex- 
periments 3, 4;] (see also below). Because QCD/ does not 
suffer from the sign problem Q , this equivalence en- 
ables us to study properties of QCD^ through the lattice 
simulation of QCD/. 



^ In this paper we consider the two-flavor QCD unless otherwise 
stated. The baryon chemical potential /ig means /ii = /i2 = 
/i = /ifl/Afc, while the isospin chemical potential /i/ stands for 
/*! = — /i2 = H = M-r/2. We also assume the degenerate quark 
mass so that the isospin symmetry is exact. 
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FIG. 1. Relations between SO(2iVc) [or Sp(2iVc)] gauge the- 
ory at finite /is (SOs or Sp^), QCD at finite fiB (QCDb), 
and QCD at finite /i/ (QCD/). A class of correlation func- 
tions in S0(2A^c) theory [U(l)fl-neutral operators; red circle] 
coincide with the counterparts in QCDb and another class of 
those (isospin singlets; blue circle) coincide with the counter- 
parts in QCD/. There is an intersection of these two classes 
in SO(2A^c), which leads to the equivalence between the coun- 
terparts in QCDb and QCD/. For the detail, see Ref. 01 • 



This equivalence has been derived by using a string- 
inspired large- iVc technique known as the orbifold equiva- 
lence fv'-^lll. (For the idea of the orbifold equivalence, see 
Sec. lIII Al ) As shown in Refs. ^, |4|], there are Z2 projec- 
tions called the orbifold projections relating SO(27Vc) and 
Sp(2A^c) gauge theoriefQ at finite (SOb and Sp^) to 
QCDb and QCD/. The relations between these theories 
are summarized in Fig. [1] The large- iVc orbifold equiva- 
lence guarantees that these theories are equivalent in the 
sense that a class of correlation functions (e.g., magni- 
tude of the chiral condensate) and the phase diagrams 
characterized by such quantities coincide, as long as the 
projection symmetry is not broken spontaneously 
This requirement is always satisfied for the equivalences 
between SOs, Sp^, and QCD/. Although the projec- 
tion symmetries relating these three theories to QCDs 
are broken spontaneously in the pion or diquark conden- 
sation phase (see Sec. Ill Bl below), the equivalence holds 
outside that region. 

The purposes of the present paper are two-fold. Firstly, 
we develop a technique of the orbifold equivalence within 
the frameworks of effective models and holographic mod- 
els which are widely used to study the properties of finite- 
density QCD. The effective models covered in this pa- 
per include the Nambu-Jona-Lasinio (NJL) model 12], 
linear sigma model (La M) [13l . Polyakov-Nambu-Jona- 
Lasinio (PNJL) model [14| - |la |. Polyakov-quark-meson 
(PQM) model nT\, chiral random matrix model (xRMM) 



[18| , and strong-coupling expansion of lattice QCD llOll . 
The holographic models include the D3/D7 model [20| 
and Sakai-Sugimoto model [2l|. We also explain impli- 
cations of the orbifold projections for lattice QCD meth- 
ods, such as the reweighting method [22,^], QCD with 
imaginary chemical potential [1, HJ] , and Taylor expan- 
sion method p5l - [27j . Then we point out that, in previ- 
ous simulation results, the phase quenching is found to 
be well satisfied at A^c = 3 for small-/z and high-T region. 
Secondly, as a consequence of the orbifold equivalence, we 
provide criteria for the validity of the phase quenching in 
a unified manner, independently of the effective models 
and lattice methods. When these criteria are satisfied, 
there is no overlap problem (see Sec. IIIIBI for details). 
Our criteria are summarized as follows]! 

1. For which quantities? — For any physical (gauge- 
invariant) observable. 

2. To what extent? — For fermionic quantities, the 
phase quenching is exact at the leading order in 
1 /Nc (planar and one fermion loop) in the large- A'c 
QCD, and in the mean- field approximation (MFA) 
in the effective models. For gluonic quantities, the 
phase quenching is exact to the next-to-leading or- 
der in 1/Nc (planar and one fermion loop) in the 
large- A'c QCD, and to the one-meson-loop correc- 
tions beyond the MFA in the effective models. 

3. In what region of the phase diagram? — The phase 
quenching is valid outside the pion condensation 
phase of the corresponding phase-quenched theory. 

Before discussing a general but rather mathematical 
derivation of the above statement based on the orb- 
ifold equivalence, let us explain a heuristic derivation 
for a specific observable at perturbative level [s^, Hlj. 
As an example, let us consider the chiral condensate 
{^tp) = (uu) + (dd). Each flavor is assumed to have the 
quark chemical potential /i„ and fj,d- At large Nc, the 
chiral condensate is dominated by one-fermion-loop pla- 
nar diagrams, with no additional fermion loops attached. 
If the flavors are not mixed, the contribution to the chi- 
ral condensate is given by the summation of each fla- 
vor, (?/''0)M„,Md = /(A^«) + /(Md) with some function /(/i). 
Here note that, as long as the ground state does not mix 
the flavors, the flavor mixing arises only through the dia- 
gram with additional fermion loop(s) which is suppressed 
at large Nc- Also note that /(/i) is an even function of /i, 
/(/i) = f{~'fi), due to the charge conjugation symmetry. 



^ The symplectic group is defined as Sp(2A''c) = {36 
\]{2Nc)\g'^ Jcg = Jc}, where Jc = —102 ® Ijv^. This is also 
denoted as USp(2Afc). 



^ We note here that, in the chiral random matrix model [28| and 
NJL model |29| . the exactness of the phase quenching was al- 
ready observed by computing the eff'ective potentials in terms of 
the chiral condensate. In this paper, we rather provide the un- 
derlying principle why this should be so, regardless of the details 
of models. As a by-product, we can also predict the exactness of 
the phase quenching in other models (LcrM etc.) which has not 
yet been pointed out, to our best knowledge. 
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Then the chiral condensate at finite hb, (■0V')Aiu=/id=M' 
turns out to be equal to the chiral condensate at finite 
M/- at the leading order of l/N^, 

= /(/^) + /(m) = + /(-Ai) = (1) 

and hence the phase quenching is exact in the large- A'c 
limit. As is found from this argument, the essence of 
the phase quenching is the flavor decoupling of chemical 
potentialso However, flavor decoupling is not satisfied if 
there is some mixing between up and down quarks in the 
ground state. This actually happens in the pion conden- 
sation phase of QCD/, as we shall see in Sec. IIIBI 

It should be remarked that the arguments based on 
the orbifold equivalence are more general. They lead to 
criteria for the validity of the phase quenching systemat- 
ically (criteria for the validity on which correlation func- 
tions and which regions of the phase diagram) based on 
the projections and the symmetry breaking patterns in a 
unified manner. Furthermore, although the proof given 
in Rcf. [4] applies to all orders in perturbation, there 
are convincing arguments that the orbifold equivalence 
holds nonperturbatively, based on the weak-coupling cal- 
culation at high density Q , effective theory analysis [s^] , 
and holographic analogue |33j . 

This paper is organized as follows. We start with re- 
viewing the sign problem and the phase quenching in 
QCDs in Secini After explaining the sign problem and 
the phase quenching in Sec. Ill A[ we show the phase dia- 
gram of the phase-quenched QCD in Sec. IIIBI In Sec. IIIII 
we argue the large- A^c equivalence which assures the ex- 
actness of the phase quenching for various observables. 
In the following sections we consider implications of this 
equivalence to effective models of QCD fSec. IIVI) . holo- 
graphic models (Sec. El), and lattice QCD (Sec. lVI[) . Sec- 
tion |VlTl is devoted to discussions and outlooks. 



II. PHASE-QUENCHED QCD 

In this section, we recapitulate the notion of the phase 
quenching. We also review the phase diagrams of the 
phase-quenched QCD (QCD/), SOb and Sp^ which are 
important for later discussions on the applicability of the 
phase quenching. 

A. Sign problem and phase quenching 

We consider mass-degenerate two-flavor QCD at a fi- 
nite baryon chemical potential fiB = Nd-i. The La- 



grangian in the Euclidean spacetime is given by 

2 

^QCD = X! ^f^i^^)^-'f + ^YM, 
/=1 

CYM^^Tr{F^.?, (2) 

where L»(/i) = 7'"D^ + m + ^7"* and D^' = + igA^'. 
ijjf is the quark field with mass m in the fundamental 
representation and — A^^T"' is the gauge field with 
T" being the SU(A^c) color generators. The partition 
function reads 

Zb = J dA [dctD{n)]^e-^™, (3) 

where S'ym is the action of the pure Yang-Mills. The path 
integral measure of the theory is positive semi-definite 
only at /is = 0. This can be understood as follows: if we 
define the eigenvalue of the operator ^^^D^^ -f /i7^ as iA„, 
it also has the eigenvalue — jA„ for A„ 7^ due to the 
chiral symmetry (i.e. this operator anticommutes with 
75). Because ^^^D^^ is anti-Hermitian, A„ are real when 
/i = 0. This is no longer true at /i 7^ where A„ are 
complex in general, A„ G C. Recalling that eigenvalues 
of D{p) appear in pairs (±zA„ -I- m) for A„ ^ 0, 

AeiD{y) = J|(iA„ -t- m)(-iA„ + m) = ^{{^1 + m^), 

n n 

(4) 

is positive semi-definite (complex) at /.t = (/x 7^ 0). The 
complex fermion determinant at /i 7^ is the notorious 
fermion sign problem, which prevents us from applying 
the Monte Carlo methods. 

One can think of the phase-quenched QCD where the 
complex phase of the fermion determinant in QCDs is 
ignored. (This is different from the usual quenched ap- 
proximation in the sense that the absolute value of the 
fermion determinant is taken into account.) This theory 
does not have the sign problem by definition, and thus, 
it can be analyzed by Monte Carlo simulations, although 
its relation to QCDs is not clear a priori. The partition 
function of the phase-quenched QCD is given bjH 

Zi = J dA\detD{fi)\^e-'^™. (5) 

The reason why we use / in the subscript is that phys- 
ically this theory corresponds to QCD at finite isospin 
chemical potential /i/ = 2fi (QCD/) @. This can be 
understood by recalling that the fermion determinant of 
QCD/ is given by 

detD(/i) • det D{-i.i) = |detD(/i)p. (6) 



Precisely speaking, the flavor decoupling is a sufficient condition 
for the exact phase quenching. Even if there is a flavor mixing as 
{i>i^)ii^,iid = /(Mu) + /(Md)+9(M«.Md) with some function 3, the 
phase quenching can be exact as long as g{lJ.u, fJ-d) = 9 (Mm ~t^d) 
is satisfied. 



^ The partition functions Zg and Zj are also denoted as Zi+i 
and Zm_it in the literature. Here 1* stands for the so-called 
conjugate quark which has the chemical potential —fi (while the 
usual quark has the chemical potential +fi) [34l . 
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Here the equality above follows from det D^—ji) = 
[det _D(/i)]*, which can be checked by using 

Jsh'^D" + m- Ai7°)7''^ = {j^Di" +m + /i7°)t. (7) 

The expectation value of an observable O in each the- 
ory is given by 

{0)b 0(detD(M))'e-^™, 

{0)i j dA C'|detL»(At)pe~^™. (8) 

Although one cannot calculate (O) b directly because of 
the sign problem, in principle one can calculate it by 
using a trivial relation 



{0)i 



(Oe 



2ie\ 



(9) 



where e'^ = det -D(/Lt) /| det -D(/i)| is the phase factor of 
the fermion determinant. This approach is called the 
phase reweighting. In practice, however, both the nu- 
merator (Oe^*^)/ and the denominator (e^*^)/ becomes 
almost zero and it is impossible to study QCD^ by using 
the reweighting method with a reasonable computational 
cost. 

Another related issue is that the phase-quenched en- 
semble (QCD/) may not have large enough overlap with 
the ensemble in the full theory (QCDs) so that the im- 
portance sampling fails; for example, it might be possible 
that the peak in the phase-quenched ensemble disappears 
because of the phase fluctuation and the tail might cor- 
respond to the peak of the full theory (Fig. [2l right). 
If this is the case, in practice the configurations around 
the true vacuum do not appear at all. This problem is 
called the overlap problem. In QCD/, as one increases 
the chemical potential, the pion condensation appears at 
some point (see Sec. IIIB|) . There QCD/ has completely 
different vacuum structure from QCD^, and the severe 
overlap problem appears. (Even outside the pion con- 
densation, the overlap problem is absent only for a class 
of observables. We will see that the orbifold equivalence 
provides us with the criteria given in Sec. HI) 

The average phase factor (e^*^)/ serves as a measure 
of the severity of the sign problem [H, H^; If the sign 
problem is mild (severe), it is close to unity (zero). Note 
however that (e^*®)/ ~ does not necessarily exclude 
(O) B = (C*)/. It happens when the phase and the observ- 
able factorize, (Oe^*^)/ = {0}i ■ (e^*^)/. If it is realized, 
we can compute some (C)b by computing (O)/ without 
suffering from the sign problem, despite a vanishingly 
small average phase factor. As we shall show below, this 
actually happens in finite-density QCD in the large-A'c 
limit and effective models in the mean-field approxima- 
tion. 




FIG. 2. Cartoons for the overlap problem. The horizontal 
axis stands for the value of an observable (e.g., the chiral 
condensate) and the blue and purple lines represent the path- 
integral weights of the full and phase-quenched theories. In 
the top figure, two distributions have large overlap. In finite- 
density QCD, it corresponds to a class of operators satisfying 
the criteria 1 in Sec.|T]in the small- region (where the pion 
does not condense in the phase-quenched theory). In the bot- 
tom figure, two configurations almost do not overlap. This 
corresponds to the large- region and/or observables which 
do not satisfy the criteria 1. 
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FIG. 3. Phase diagram of phase-quenched QCD at finite /i 
(or QCD at finite = 2/i). Figure taken from Ref. [3]. 



B. Phase diagrams of QCD/, SOb, and SpB 

In this subsection we discuss the phase diagrams of 
the phase-quenched QCD (QCD/), SOs and Sps, with 
particular emphasis on the symmetry breaking patterns. 
The phase diagram of QCD/ in the T-^i plane was first 
investigated in Ref. 6]. It was shown recently Q that, in 
the large- iVc limit, SOs and Sp^ have exactly the same 
phase structures as QCD/0 Even at finite iVc, they are 
qualitatively the same. The phase diagrams of QCD/ 



® In the large- iVc limit and the chiral limit, the phase diagram of 
QCD/ also coincides with that of QCD at a finite chiral chemical 
potential /is |37| . where /is is defined as the chemical potential 
associated with the U(l) axial charge. 
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FIG. 4. Phase diagram of SO(27Vc) gauge theory at finite /i. 
Figure taken from Ref. 0|. 

and SOb are depicted in Fig. |3] and Fig. 21 respectively. 

We first consider the ground state of QCD/ at T = 0. 
As Hi is increased, what happens first is the excitation of 
the Hghtest particle with the isospin number, the charged 
pion. At /J,/ > m,r (m > "Itt/S), where ttItt is the pion 
mass in the QCD vacuum, the excitation energy of a 
charged pion, — ^J■I, becomes negative and it is ener- 
getically favorable to excite it. Since a pion is a boson, 
this is the Bose-Einstein condensation (BEC) of pions 
which breaks \J{1)l+r down to Z2. 

In the opposite limit, at sufficiently large /i/, inter- 
actions around the Fermi surface become weak because 
of the asymptotic freedom, and fundamental degrees 
of freedom are quarks and gluons. The one-gluon ex- 
change interaction between quarks near the Fermi sur- 
face leads to an attractive interaction in the color singlet 
channel. According to the Bardeen-Cooper-Schrieffer 
(BCS) mechanism, the Cooper pairing (d'^^u) is formed 
[ijlll Because the quantum numbers of condensates and 
symmetry breaking patterns are same in the BEC and 
BCS regimes, it is plausible that these two regimes are 
smoothly connected without any phase transition. This 
is the so-called BEC-BCS crossover. 

As T is increased, the pion condensation is melted away 
and \J{1)l+r symmetry is restored. In the BEC regime, 
the critical chemical potential fic for the pion conden- 
sation is not so sensitive to T at low temperature. In 
the BCS regime, the critical temperature Tc should be 
an increasing function of fi, since the phase space near 
the Fermi surface for the pairing becomes larger with in- 
creasing III. A natural scenario to continuously connect 
these two regimes is shown in Fig. |31 as first proposed in 
Ref. 0. 

Next let us consider SO^ . The Lagrangian of the the- 
ory is given by Eq. ([5]) with the color generators re- 
placed by those of the gauge group S0(2iVc). A crucial 



Precisely speaking, the one-gluon exchange interaction does not 
distinguish between the condensates {d'ysu) and (du). The con- 
densate {d75u} is favored by the nonperturbative instanton- 
induced interaction 



difference from QCD is that there is no distinction be- 
tween fermions in the fundamental and antifundamental 
representations (in the vacuum) because the gauge group 
is real. For this reason, mesons in this theory are not 
necessarily neutral under U(1)b; baryon- number charged 
mesons (baryonic mesons or antibaryonic mesons) out of 
two quarks or two antiquarks also arise. 

In order to identify the lightest baryonic meson, let 
us consider the symmetry breaking pattern of the the- 
ory. When m = /iB = 0, the chiral symmetry, which 
at first sight looks SU(iV/)i x SU(iV/)i? x U(1)b, is en- 
hanced to SU(2A^/), where Nf is the number of flavors. 
This enhancement of chiral symmetry originates from the 
fact that there is no distinction between left- and right- 
handedness. One can actually write the fermionic part 
of the Lagrangian manifestly invariant under SU(2iV^) 
using the new variable ^' = (V'l, <^24'*r)^ '■ 

Ci = i-^^a^D^-^, (10) 

where = (tJfc,— 112) with the Pauli matrices 
{k = 1,2,3). The SU(2A^/) chiral symmetry is spon- 
taneously broken to S0(2iV/) by the formation of the 
chiral condensate {'ip^')^ which gives rise to 2N'j + 
Nf — 1 Nambu-Goldstone bosons parametrized by the 
coset space SU(2A^/)/SO(2A^/): neutral mesons Ua ~ 
ipJbPa^^ baryonic mesons (or simply diquark) = 
4''^C'j5Qs4' ^-nd antibaryonic mesons (or antidiquark) 
~ "'P^C^lzQsiJ'* 1 where Pa are traceless and Hermi- 
tian Nf x Nf matrices. Pa = (a = 1, 2, • • • , - 1), 

and Qs are symmetric Nf x Nf matrices, — Qs 
(5 = 1, 2, • • • , Nf{Nf + l)/2), in the flavor space. In the 
vacuum, their masses are degenerate due to the unbro- 
ken flavor symmetry, mn = = m^t- (At large iVc, 
because of the orbifold equivalence, their masses also co- 
incide with the pion mass ttItt in QCD.) This degeneracy 
is resolved if we turn on nonzero /i which explicitly breaks 
the SO(27V/) symmetry down to SU(iV/)y x U(1)b. 

At r = 0, as we turn on /x^, the lightest particle 
with the baryon number, the diquark E5, condenses for 
/i > mY,/2. This is the BEC of diquarks where U(1)b 
symmetry is broken to Z2. At sufficiently large /x^, on 
the other hand, the one-gluon exchange interaction be- 
tween quarks near the Fermi surface is attractive in the 
color symmetric channel and gives rise to the pairing of 
the form {tp'^Cj^Qsip) 0. This BCS pairing breaks the 
same symmetry as the BEC at small /is, and it is again 
natural to expect the BEC-BCS crossover for fj, > ms/2 
along the fi axis. 

In the same manner, one can also obtain the phase 
diagram of Sp^. (For further details, see Ref. Q-) 

III. PHASE QUENCHING IN LARGE-A/^c QCD 
A. Orbifold equivalence 

In this subsection, we briefly review the large- A''c orb- 
ifold equivalence [THllj and apply it to QCD and QCD- 



6 



like theories [13111. Thereby we estabhsh the exactness 
of the phase quenching in the large- A^^c limit (for details, 
see Ref. The relations between QCD and QCD-like 
theories through the orbifold projections are summarized 
in Fig. [H 

The idea of the orbifold equivalence is the following: 
first we choose the discrete symmetry P (subgroup of 
gauge, flavor, or spacetime symmetry) of the original 
theory called the parent. We then throw away all the 
degrees of freedom not invariant under P. This proce- 
dure is called the orbifold projection. After the projec- 
tion, we obtain a new theory called the daughter. The 
orbifold equivalence states that, in the large- TVc 't Hooft 
limit where the 't Hooft coupling g^Nc is kept finite, cor- 
relation functions of operators O^^^ {Ai^,ip) invariant un- 
der P in the parent (called neutral operators) agree with 
those of the operators ©('^^(AP^j, ^/-P^j) that consist of 
projected fields in the daughter: 

{0^f^0iP^---), = {0['^0i''^ ■■■),. (11) 

Here coupling constants should be appropriately related; 
for example, for the equivalence between QCD^ with 
SU(7Vc) gauge group and SOs with S0(2A^c) gauge 
group, which we shall consider below, we take 

5iu = 9io- (12) 

The field theoretic proof to all orders in the perturbation 
theory was given by Bershadsky and Johansen [l^ and 
nonperturbative proof in certain gauge theories was given 
by Kovtun, Unsal, and Yaffe [lli]. For QCDb, QCD/, 
SOs, and Sps, a couple of evidences of nonperturba- 
tive equivalence were also provided by the weak-coupling 
analysis at high density limit |^] , chiral perturbation the- 
ories ^11, chiral random matrix models [1], and holo- 
graphic models [33l|. 

Here let us take SOb as a parent and consider the 
projections to QCD^ and QCD/ independently^ We 
identify the Z4 discrete symmetries of SO b generated by 
= -ia2 (g) Itv, G S0{2Nc) and w = e*''/^ g u(l)/3, 
where 1 jv is an x identity matrix^ We require the 
gauge field ^^^^ and the fermion to be invariant 
under the following Z2 transformation embedded in the 
gauge and U(1)b transformations @, 

^^^b — iJc)aa'A^^a'b'i-Jc ^)b'b, (13) 



For an earlier work of the orbifold projection from SO(2Afc) to 
SV{Ni:) gauge theories, see Ref. 1381 . See also Ref. [3911 where 
phase diagrams of QCD-like theories with different matter con- 
tents at small have been studied in the context of the orbifold 
equivalence. 

Here Jc is chosen to satisfy the regularity condition, Tr(J") = 
when J" 7^ itl2]Vc . This condition is necessary for the derivation 
of the orbifold equivalence [l^ . 



Under these projection conditions, QCD^ is obtained as 
the daughter. In order to see it, we decompose the gauge 
field and fermion field of the parent SOb as 

(15) 

where the gauge fields with the superscript A (S) are 
Nc X Nc anti-symmetric (symmetric) matrices and ^ and 
C are A^c-component fermions. Under the Z2 symmetry, 
and D^, are even while and are odd, and the 
orbifold projection sets B^ = = 0; we can also see ^ 
and C are even and odd under Z2, respectively. Therefore, 
the daughter fields after the projection are 

^r=i^s-fl), ^r = (|) (16) 

After a unitary transformation using the matrix 
it can be rewritten as 

p /IProj p-l — ^("^m)^ ^ ^ p,/,proj _ \ 

P\ P '[ (AU)J' -[rj' 

(18) 

where A^ = + iA^ is a U(Afc) gauge field and 
-0^ = "v/Z^. Since the difference between U(A"c) and 
SU(A''c) is a l/A"^ correction and is negligible at large Nc, 
the daughter theory can be regarded as QCD^ given by 
Eq. From this orbifold projection, we have the equiv- 
alence between SOs and QCD^. However, the U(1)b 
symmetry, whose Z4 subgroup is used for the projection 
of the fermion in Eq. (|14p . is spontaneously broken to 
Z2 in the diquark condensation phase (the BEC/BCS re- 
gion in Fig. in ; the equivalence is valid only outside that 
region. 

One can also construct the projection from SOs to 
QCD/ by choosing another Z2 symmetry [S, [4], 

Al'^ab={Jc)aa'Al%{J-'),.,, (19) 
i'l%-^{Jc)aa'i^l%r[Jr')i'f^ (20) 

where Ji = —icf2 ® generates Z4 subgroup of SU(2) 

isospin symmetry and the projection condition for the 
gauge field is the same as Eq. To see how QCD/ 

can be obtained through the projection, we decompose 
the flavor 2A^^-component fundamental fermion into two 
A^y-component fields, 

V^^o = ^,), (21) 

with i and j being the isospin indices, and furthermore 
decompose 2A''c color components to two sets of A^c com- 
ponents, 

*=(<■)■ *=(<:)■ 
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If we define ^± = {^±iC)/V2, ip± = (V'i T«V'i)/%/2 and 

X± = (V'i i ) /v^, the fermions ip± survive but x± 
disappear after the projection ([20]) . Because (p± couple 
to [A^^)'-^ and A^^ respectively, the Lagrangian of the 
daughter theory is now 



C 



QCD, 



Vi/±(7''^/. + m±M7')V'^±, (23) 



E 



where ■0+'^ = V^'^- and ■0-^ = "\/2'^+- This theory is 
QCD/. In this case, the isospin symmetry used for the 
projection of the fermion is unbroken everywhere, and so 
the orbifold equivalence holds including the BEC/BCS 
region of the phase diagram. Therefore, through the 
equivalence with SOb , we obtain the equivalence between 
QCDs and QCD/ outside the BEC/BCS region of the 
latter; the phase quenching is exact for neutral sectors in 
this region. The same conclusion can be reached through 
the equivalence with Sp^ [2] . 

A few remarks are in order here. Firstly, not all the op- 
erators coincide. In the parent theory, only the operators 
invariant under the projection symmetry P are related to 
the counterparts in the daughter. In the daughter, not all 
the operators are obtained from the parent through the 
projections. As an example, consider the fate of neutral 
pious and (anti)diquarks of the parent SOs after the orb- 
ifold projections (see Tab. |T|. The projection to QCD^ 
maps U(l)B-neutral pious of SOb to pious in QCD^, 
TT*^, 7r+ and 7r~, and throws away (anti)diquarks. On the 
other hand, the projection to QCD/ sends (isospin singlet 
part of) U(l)B-neutral pions, diquarks, antidiquarks, to 
7r°, 7r+ and tt", respectively. Therefore the diquarks in 
S0(2A^c) theory and 7r+ in QCD/ have the same mass 
TOjr in the vacuum and the same excitation energy at any 
jjL (at T = 0). In the same way as 7r+ condenses in QCD/ 
at = m7r/2, the diquark condenses at = TOtt/S (see 
Fig. [3] and Fig. 0]). Note that the charged pions i:^ in 
QCDb and QCD/ have different origins in SO^ and they 
do not correspond each other. 

Secondly, note that two projections and (pOt are 
equivalent when /i = as they should be. Both are a Z4 
subgroup of the flavor symmetry which mixes two Ma- 
jorana flavors. Once {^b or /^/) is turned on, they 
are not equivalent. The flavor symmetry Ji used in (j20l) 
is essentially the same as Jc, and the proof in Ref. [lOj 
can be repeated straightforwardly [4] ; the only difference 



TABLE I. Some examples of the correspondence between 
SOs, QCDj and QCDb through the orbifold projections. See 
the text for detail. 





Order parameter 


Elementary excitations 


SOs 
QCD/ 
QCDs 


(l/'-i/'>SOB 
(V'^>QCD^ 
{•(/'^>QCDb 


Neutral pions 

^» 
+ - 


Diquarks 

X 


Antidiquarks 
7r~ 

X 



is some color-index loops which are replaced by flavor- 
index loops. On the other hand, Z4 € U(l) b used in 
Eq. (|I4p is different and the proof in Ref. holds only 
for planar diagrams with at most one fermion loop. Be- 
cause fermion loops are suppressed by the factor Nf/Nc, 
the equivalence through the projection (fT4|) to QCDs 
holds in the 't Hooft large-iVc limit (Nc — >■ 00 with Nf 
fixed) while the one through the projection to QCD/ pp)) 
holds also in the Veneziano large- A^c limit {Nc — > 00 with 
Nf/Nc fixed) 

The above second remark has an implication for the 
1/A^c corrections Compare QCDs and QCD/. In 
the 't Hooft large-A^c limit, expectation values of glu- 
onic operators trivially agree because the fermions are 
not dynamical. Now consider finite- A'^ci say Nc = 3 and 
Nf = 2. Then the largest correction to the 't Hooft limit 
comes from one-fermion-loop planar diagrams, which, as 
we have seen, do not distinguish /is and /i/. Therefore 
the difference of expectation values of gluonic operators 
is at most (Nf/Nc)^ (two- fermion- loop planar diagrams). 
In particular, the deconfinement temperatures, which are 
determined by the Polyakov loop, agree up to corrections 
of this order. A similar observation was made in Ref. (3l| 
by a perturbative argument. 

Note that the 1 /Nc correction can become larger in the 
confining phase, because of thermal excitations of pions, 
resonances, and baryons, which large-A^c arguments do 
not take into account; baryon gas in QCD^ is quite dif- 
ferent from pion gas in QCD/ [36|. On the other hand, 
for T > Tc, fundamental degrees of freedom are decon- 
fined quarks and gluons rather than baryons and mesons, 
where the difference between QCDs and QCD/ becomes 
much smaller and the large- A'c equivalence may be well 
satisfied even at Nc — 3. It can indeed be confirmed 
numerically, as we will see in Sec. IVIl The results there 
show that the phase quenching is a very useful tool for 
the study of the chiral transition. 



B. Implications for the phase reweighting 

In the phase reweighting method, one calculates ob- 
servables (e.g., the chiral condensate) by using the QCD/ 
ensemble and by taking into account the phase factor. 
There, the absence of a severe overlap problem is implic- 
itly assumed. But how can one justify this assumption? 
Actually we can obtain a very simple answer by com- 
bining the orbifold equivalence and a simple symmetry 
argument: any observable coincide to 0{Nf /Nc) imply- 
ing that there is no overlap problem for these quantities 
to this order. 

At first it may sounds implausible; for example, the 
baryon number density takes nonzero value in QCD^, 
while it is zero in QCD/, but then how the phase quench- 
ing can work? 

The point is that, although the partition functions of 
QCD/ and the phase-quenched QCD are equivalent, ob- 
servables are not necessarily the same. For example. 
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in QCD/, the propagators of up and down quarks are 
_D^^(+/i) and D~^{—^), while in the phase-quenched 
QCD, both of them are D^^{+ii). (In the terminology of 
the lattice simulation, we generate the configuration by 



using QCD/, but the measurement is done by using the 
same operators as QCD^.) As a result, the expressions 
for the chiral condensate and the baryon/isospin density 
in each theory are given by 





qcdb 


QCD/ 


Phase quenched QCD 




2(TrZ?-i(M)>B 


(TrL»-i(/Li) + TrZ?-i(-/i))/ 


2{^TD-\^i))I 




2(Tr 707^-1 (/x))b 







{ni) 





(Tr7"L>-i(^) - Tr7''D-i (-//)) 7 






Because {TtD-'^{p))i = (Tr (-//))/ due to the 
charge conjugation invariance, the chiral condensate in 
QCD/ and that in the phase-quenched QCD take the 
same value. Therefore, the orbifold equivalence states 
that it remains unchanged by the phase quenching. For 
the baryon density, we first note that (n/)/ — {tib) i be- 
cause (Tr7°i:)-i(-^))/ = -(Tr7°D-i(^))/. By com- 
bining it with the orbifold equivalence, {nB)B = 
we conclude that the phase quenching is valid for the 
baryon density. Note also that the isospin density is 
trivially zero in QCD^ and in the phase-quenched QCD 
it is {TT-f°D-'^{^i) - TT-f°D-'^{^i))i = 0, and the phase 
quenching is valid. The same argument is applicable to 
other fermionic observables too. 

The exactness of the phase quenching in the large- A^'c 
limit can be rephrased that the overlap problem is 1 /Nc- 
suppressed. Note however that the phase reweighting is 
doable only at small volume and Nc] the average phase 
factor is exponentially suppressed at large volume and/or 



C. Relation to the quenched approximation 

Let us consider the quenched approximation in lat- 
tice QCD, in which the fermionic observables are cal- 
culated by using the gauge configurations generated in 
the quenched QCD (i.e., the pure Yang-Mills theory). 
This approximation is believed to become exact in the 't 
Hooft large- iVc limit where dynamical fermion loops are 
suppressed. Does this approximation make sense at finite 
A*? 

For concreteness, let us consider the chiral condensate. 
It is calculated as 



('/^/V'/)YM = (Tr(^7'(AA^/) 



YM 



(24) 



where the propagator DJ^{A, ^) is a function of the 
gauge field and the path integral is taken by using 
the pure Yang-Mills action. As long as one considers a 
perturbation around the trivial vacuum, it contains all 
one-fermion-loop diagrams and gives the right answer at 
large Nc. However whether it is correct at nonperturba- 
tive level is nontrivial, and actually it fails at finite fi - 



although QCDb and QCD/ are identical in the quenched 
approximation as we have seen in Sec. HI QCDg and 
QCD/ actually have completely different phase diagrams. 
Then what is wrong with the quenched approximation? 

The expectation values in QCD^ and QCD/ are writ- 
ten in terms of the quenched QCD as 



{'^f^f)BJ = 

Tr(D-\A,ijLf)) ■Y[^AelDf{A,^Xf] 



YM 



n,detD/(A/^/) 



(25) 



YM 



They agree with the value under the quenched approxi- 
mation when the following factorization holds. 



l:,(D-\A,^is)) -ndetDKA 



YM 



Tr 



{D-\A,^Jif) 



YM 



\[dciDf{A,^lf) 



YM 



(26) 



This factorization should be distinguished from the usual 
one which follows from the 't Hooft counting; although 
a finite number of traces factorize, the factorization is 
not valid for determinants. It is plausible that the fac- 
torization takes place in QCD/, because the quenched 
approximation exhibits the pion condensation, as is ex- 
plicitly demonstrated within the chiral random matrix 
model [s^l- However, the factorization obviously fails in 
QCDs for > m7r/2 at T = where the pion conden- 
sation should not occur. Although we have not proven 
in this paper, the equivalence between QCD^j and the 
quenched QCD outside the pion condensation region at 
large would be useful because the quenched QCD is 
numerically cheaper. It is also interesting to calculate 
the chiral condensate and the baryon/isospin density by 
using lattice configurations with dynamical fermions at 
= 0. 



IV. PHASE QUENCHING IN EFFECTIVE 
MODELS OF QCD 

As we have seen, the phase quenching is exaci 
the large- TVc hmit for a class of observables outside 
pion condensation phase of the phase-quenched the 
(QCD/). Therefore effective models of QCD, if desc: 
underlying QCD properly, should exhibit the same la: 
Nc equivalence. Then the phase quenching is expec 
to be exact in the mean-field approximation (MFA 
the models, since MFA corresponds to the leading-O] 
in the expansion of QCD. 

In this section, we show this statement by explic 
constructing orbifold projections within several mo^ 
frequently used to study the phase diagram of Ql 
including the Nambu-Jona-Lasinio (NJL) model, 
ear sigma model (LcrM), Polyakov-Nambu-Jona-Las 
(PNJL) model, Polyakov-quark-meson (PQM) mo 
chiral random matrix model (^RMM), and strc 
coupling expansion of lattice QCD. (The orbifold ] 
jections of the xRMM were already given in Ref. [1] , 
we include it here for completeness.) We then anal 
cally demonstrate the exactness of the phase quencl 
in NJL model Hi and xRMM [I,!!!], for which the 
fective potentials are known. 



A. Mean-field approximation 

In order to apply the large- iVc equivalence to the 
fective models, let us set up the 1/iVc-counting sch( 
in the models, so that the right powers of 1/Nc in Q 
are reproduced. As an example, we consider the I 
model (see Sec. IIVBI for the Lagrangian). The qr 
field ip has Nc colors so that a closed color loop g 
a factor of Nc. (Here Nc is treated as a variable and 
will be taken to A'^c = 3 at the end of calculations.) We 
need to make the counting scheme such that each fla- 
vor loop gives a suppression factor of l/Nc- Then the 
coupling constant of the four-fermi interaction should be 
taken as 0(A~^), and furthermore, the form of possi- 
ble four-fermi interactions are restricted: the interaction 
of the form {^af'4'ag){'4'bg'4'bf) showu in Fig. [SJa) is al- 
lowed, but {■tpafipaf)itpbg'4'bg) in Fig. [Sfb) is not, where 

b (/, g) are color (flavor) indices. This is because the 
one-flavor-loop diagram in Fig.lHKb) derived from the lat- 
ter interaction is not suppressed in compared with 
the diagram with no flavor loop in Fig. [S] (al) derived 
from the former. Once we exclude the interaction of the 
form {■ipafipaf)i^bg'ipbg), the right 1 /Ac-counting follows 
and we can use the same proof of the orbifold equivalence 
in the NJL model as the large- A^ QCD in Ref. 

Under these conditions, we now recall the strategy 
to compute the effective potential of the NJL model. 
We first perform the Hubbard-Stratonovich transforma- 
tion by introducing auxiliary fields corresponding to 
the fermion bilinears, a a = {G/Nc)tpTAip and tt^ — 
{G/Nc)ipi^5TAip, where ta are the U(2) fiavor generators. 



a a 

b b 
(a) 



af 



9^ 



> < 



af 



gb 



(b) 



FIG. 5. Two types of interactions: (a) {il)af'ipag){ipbg4^bf) 
originating from the one-gluon exchange interaction and (b) 
{tf'af'4'af){'4'bg'4'bg) Originating from the instanton-induced in- 
teraction, with a, b (/, g) being color (flavor) indices. The 
solid and dotted lines denote color and flavor lines, respec- 
tively. 




(a1) 



(a2) 




FIG. 6. Loop diagrams induced by the two types of interac- 
tions (a) and (b) in Fig. [S] (al) 0{N~^) with no flavor loop, 
(a2) OiN-'^) with one flavor loop, and (b) 0{N^ ^) with one 
flavor loop. 



We then integrate out fermions to obtain the partition 
function 



(27) 



Here /(cr, tt) is the bosonized action 



/((TA,7r^) = Nc 



Trlogi?+i 



(28) 



with D = 7^9^ -1- 2{aA + tta)- This effective action de- 
scribes a theory for bosonic fields (mesons) a a and tta- 
Because there is an overall factor Nc in the action, the 
expansion of Z [or /((7^,7r^)] in terms of 1/Nc is equiv- 
alent to the expansion in terms of meson loops In 
particular, the leading order in 1/Ac corresponds to the 
saddle-point approximation, or the MFA where the auxil- 
iary fields are replaced by the expectation values (i.e., the 
mean fields) In the language of many-body physics. 



A similar discussion on the 1/N expansion in the context of con- 
densed matter physics, e.g., fermions at unitarity, can be found 
in Refs. [JT, 42| where A'^ is the number of species of fermions. 
In Rof. 41] 1/N corrections are found to be numerically small 
by Monte Carlo simulations. 
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the MFA corresponds to the Hartree approximation for 
the quark self-energy and to the random phase approxi- 
mation for the four-fermi interaction (see, e.g., Ref. p3l|). 
In order to go beyond the MFA, we have to take into ac- 
count n-meson-loops (n g N) order by order which give 
l/N'^ corrections to the MFA [4i |. 

Since the LcrM, PNJL model, and PQM model can 
be regarded as simple generalizations of the NJL model, 
as we will argue below, we can develop similar counting 
schemes to apply the orbifold equivalence to them. In the 
case of the xRMM, we can use the large- equivalence 
by identifying the the size of the matrix A^ as a variable 
which is finally taken to be infinity corresponding to 
the thermodynamic limit. 

In the following, we will argue the orbifold equiva- 
lence in each model. It then predicts the exactness of 
the phase quenching in the MFA outside the pion con- 
densation phase of the phase-quenched model (which we 
denote "model/"). 

B. Nambu— Jona-Lasinio model 

We first consider the Nambu- Jona-Lasinio (NJL) 
model [l3 which captures the physics of chiral symme- 
try breaking in QCD (for reviews, see Refs. [1^, Ey]). In 
order to simplify the discussion, we consider the chiral 
limit TO = so that we maintain the full chiral symmetry 
of the model. The generalizations to include nonzero m 
is straightforward. 

The starting point is the Lagrangian with the U(A'c) 
color current interaction with Nf flavors, 

/:njl = V^/ {rO, + ^/ - ^/^2j'^^2^ (29) 

where = V'/T/.^uV'/ and are the U(A^c) color 

generators and summation is taken over repeated indices. 
The coupling constant G is taken to be of order N^. One 
rewrites it keeping only the interactions in the scalar and 
pseudoscalar channels after Fierz transformations: 

-Cnjl = i' {l^df, + fif-f^) ^ + Lint, (30) 

G - - 
Ant = [(^/V'/' )('/'/"/'/) + ('/'/i7^V'/')(V'/'«7^V'/)] • 

(31) 

In the Lagrangians ([29]) and pOI) . the invariance under 
U(A^c) gauge symmetry and U(A^/)i, x \J{Nf)ji flavor 
symmetry are manifest. Here we ignore the effect of in- 
stantons or the \J{1)a anomaly which explicitly breaks 
the U(1)a symmetry, because it is subleading in l/A^cE3 



To understand this statement, let us compare the couplings of 
the one-gluon exchange interaction Gqge a^nd the instanton- 
induced interaction Gjusti where Gqge E'nd Ginst a^r^ defined as 
the coefficients of multi-fermi interactions [see Eq. H41|l for the 



from the viewpoint of the orbifold equivalence, there is 
no reason for the exactness of the phase quenching if we 
take into account the l/A'c-suppressed instanton effects. 
However, as we shall see in Sec. IIVBI even if we incorpo- 
rate them, the phase quenching for the chiral condensate 
turns out to be exact within the NJL model at the level 
of MFA. 

For S0(2A'c) theory, we can construct the correspond- 
ing NJL model in the same manner, by starting with 

where 

- (^/L, ^^fl) = {^fL,Gi^fRf) , (34) 

(F = 1, 2, • • • , 2A^/) are 2A''c component fermions and the 
current J^^'-'^ is defined by 

Jj^r = ^n.Tio^f = ^F^.Tio^F ■ (35) 

The invariance under S0(2A''c) gauge transformation and 
SU(2A^/) flavor rotation is manifest at this level. After 
the Fierz transformations and concentrating on the in- 
teractions in the scalar and pseudoscalar channels, La- 
grangian reduces to 

G - - 

(36) 

The baryon chemical potential in S0(2A'c) theory cor- 
responds to /ii? = {^fL,—lJi-fFi) = (+Mi~/^)- Note the 
minus sign in front of i^fR, which arises because of the 
charge conjugation. Due to this sign, the chiral symme- 
try is exphcitly broken to S'\]{Nf)L x SU(Af/)i?. 

From the S0(2A''c) NJL model, it is possible to get 
SU(Afc) NJL model at finite ^lB (NJL^) and SU(Afc) NJL 
model at finite /x/ (NJL/). The projections are the same 
as those used for the fermion for S0(2A''c) gauge theory 
to QCD/3 or QCD/, (fTil) and (PU]): 

,Ar=U;(J.)aa'^fP, 



form of the instanton-induced interaction]. One then finds that 
GoGE ~ ~ in the 't Hooft limit (where g is the QCD 

coupling constant), as is consistent with the 1/Nc counting for 

G/Nc ~ in Eq. ((29J. On the other hand, Ginst ~ N~^^ , 

since Ginst is related to the r]' meson mass as 

m^, Gi^^t-:^ ^ N-\ (32) 
J 

1/2 

where a ^ Nc is the chiral condensate, /^/ ^ Nc is the decay 
constant of r]' , and the relation m^, ~ follows from the 

Witten-Veneziano formula. Therefore, the instanton-induced in- 
teraction is suppressed compared with the one-gluon exchange 
interaction for N-f > 2 in the l/Nc counting. 
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respectively. Because the fermion kinetic term of the 
NJL model is the same as that of SO(27Vc) gauge theory 
except that the gauge field is now absent, the projection 
conditions above lead to the fermion kinetic term with 
/is or fij (without the gauge field), respectively. 

The projection of the four-fcrmi interaction is also sim- 
ple. Since the the current interaction of S0(2iVc) NJL 
model is mapped to the current interaction of SU(A^c) 
NJL model, if we concentrate on scalar and pseudoscalar 
sectors and eliminate all of the others, the resultant in- 
teractions in Eqs. ([55]) and ([50]) must correspond to each 
other. From these orbifold projections, NJL^ and NJL/ 
are equivalent in the MFA outside the pion condensation 
phase of the latter. 

In order to check the exactness of the phase quench- 
ing, let us look at the free energy of two-fiavor NJL model 
with quark chemical potentials fj,u and fid- This was al- 
ready computed in Ref. [2^. In the absence of the pion 
condensation, the free energy in the MFA is given by 



±,/=u,d 



(37) 



where Ej = ^p^ + Mj and Mf = to/ - {AG/Nc)(Tf. 

From the expression above, the effective potential satis- 
fies the relation 



i^NJL(M-B)U7=0 — f^NJL(M/)Uf3=0- 



(38) 



Thus, the chiral condensates c/ obtained from the gap 
equation 



NJL 



da* 



= 0, 



(/ = u,d) 



coincide, 



O-/(mb)U7=0 = Cr/(/X/)|^g=o- 



(39) 



(40) 



Therefore, the phase quenching for the free energy and 
chiral condensate is exact in the MFA. 



Instanton-induced interaction 

The instanton-induced interaction, also known as the 
Kobayashi-Maskawa-'t Hooft interaction ^|, is used in 
the practical calculations of the NJL model at the level of 
the MFA (for the treatment of the instanton-induced in- 
teraction in the MFA, see, e.g., Ref. [45l|). The instanton- 
induced interaction is suppressed compared with the one- 
gluon exchange interaction in the 1 /Nc counting (see the 
footnote in Sec. IIV Al) . However, even if we include this 
interaction, the phase quenching is still exact. This is 
because contributions of different chemical potentials, fiu 
and fid, decouple in the MFA, independently of the forms 



of interactions, as we shall show explicitly below [see 
Eq. O]. 

The instanton-induced interaction has the form 

-Cinst = -Ginst det [^f{l + 75)?Ag] + h.c, (41) 

f,g 

which respects SU(iV/)i x SU(iV/)fl x U(1)b, but breaks 
V{1)a explicitly. For two flavors, this can also be rewrit- 
ten as 



+ (V5*75rV)2], 



(42) 



where t° are the SU(2) flavor generators. The free energy 
in the MFA is computed as [29[ 

^^NJL+inst(M«j Md, r) = ^^NJl(Mm, Md, 7") + "^G inst'^u'^ d , 

(43) 

where fljqjhifJ'u, fJ-djT) is given by (l37t with Ef replaced 

by Ef =. 



ip^ + Mj, and Mf ^ mf - {4G/Nc)af - 

2Gi„st(T/- with /' ^ /. 

Despite the presence of mixing terms, such as ~ cr„(7rf, 
one flnds the potential still satisfies the property 
([55)1 . This can be understood as follows: the effective 
potential in the MFA is given by a summation of ring 
diagrams where a number of chiral condensates are at- 
tached to one central fermion loop. It is only this fermion 
loop which has the chemical potential dependence, or 
fid- The contributions of /z„ and fid are decoupled, which 
leads to the property ([55)1 . 

Therefore, the phase quenching for the free energy and 
chiral condensate, Eqs. ([551) and (HU)) . is exact in this case 
as well, as pointed out in Ref. [29,]. 



C. Linear sigma model 

The linear sigma model (LuM) , also known as the Gell- 
Mann-Levy model [13], is another model that describes 
chiral dynamics of QCD. Essentially, this is a bosonized 
theory of the NJL model with adding potential terms for 
meson fields. The Lagrangian of its SU(iVc) x \J{Nf)L x 
\J{Nf)f(_ symmetric generalization is given by 



C 



LctM 



(44) 



where 



^B = ^ [{d^T^Af + {d^cTAf] + Ui<JA, tta), (45) 



^-,(46) 



and 



U{aA,7TA) = [<JA^ + 7rA^ - (iVc/)']' - i?ao,(47) 
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where ga ~ '4'taiP (ctq ~ V'V') ^.nd tta ~ ipij5TA'4' with 
the U(-/Vj^) flavor generators ta- Note that we have in- 
cluded flavor nonsinglet scalars and flavor singlet pseu- 
doscalar to the conventional LtrM (see, e.g., Ref. (49| ) 
to maintain the U(iV/)L x \J{Nf)ji symmetry. The pa- 
rameters g,\,f and H are taken to be of order iV°F^ 
With this normalization, (ctq) ~ in the chiral symme- 
try broken phase. Also the coupling constants involving 
a A and tta, which are physically identifled with fermion 
bilinears (mesons), correctly reproduce the usual power 
counting in the large- A'^ QCD. 

To argue the orbifold equivalence, we consider the 
S0(2iVc) gauge group counterpart of the LctM, whose 
Lagrangian is given by 



£SO 



LcrM 



SO 



(48) 



where 



so 



N, 



-U{(TA,i^A,d\,dj^), 



(49) 



+ ^4''^C{TAd^ + i-f^TAd+)4> 
+ j^^CiTAd^ + n'^TAd+)^^, 



and 



U{a,7f) 



I4P 



\d 



A? 



{Ncf? 



(50) 



(51) 



where d\ 



^"^Cij^TAip and ~ tp'^CrAip (indices ± 
denote the parity) and is defined in Eq. (p4| . 

The orbifold projection from S0(2iVc) LctMb to 
SU(iVc) LaMs can be defined as follows: the projection 
for fermions is the same as Eq. (IT4l) . 

^fo = c.(J,),,,7^fp^ (52) 



break the projection symmetry (i.e., outside the diquark 
condensation phase). 

To understand the orbifold equivalence in the bosonic 
sector, we consider a neutral meson-meson scattering be- 
tween S0(2iVc) LctMb and SU(A^c) LcrMs as an exam- 
ple. (A similar argument can be found in Ref. [s^].) The 
generalizations to general scattering amplitudes and to 
the case between S0(2A^c) LaMj and SU(Afc) LctMb are 
straightforward. 

First note that a meson loop is absent in the large- A'c 
limit (MFA); we can concentrate on tree- level scatterings 
where the external lines are neutral mesons. The external 
legs are neutral pions in S0(2A^c) LctMb and are 7r°s in 
SU(Ac) LuMb (see Tab. H]). Because the neutral-meson 
coupling constants are taken to be the same between 
the two, a possible difference of the scattering amplitude 
comes from appearance of charged mesons (diquarks and 
antidiquarks) in the internal lines of the scattering dia- 
gram in S0(2A^c) LctMs, whose counterparts do not exist 
in SU(A^c) LctMb- However, this is impossible due to the 
conservation of the global Z2 charge, as long as the Z2 
symmetry is not broken spontaneously: when the exter- 
nal legs are neutral mesons, mesons must also be neu- 
tral in the internal lines of the diagram. One might still 
suspect that a pair of diquarks (antidiquarks), which is 
neutral under Z2 , could appear in the diagram. But this 
necessitates a meson loop and is suppressed in the large- 
Ac limit. Hence, the neutral meson scattering amplitude 
must agree between the two. 

From the equivalence in both fermionic and bosonic 
sectors, the equivalence holds in the full theories be- 
tween S0{2Nc) and SU(A'c) LaMs in the MFA, as long 
as the projection symmetry is unbroken. One can simi- 
larly show the equivalence between S0(2A"c) LctMb and 
SU(A^c) LcrM/. Therefore, the phase quenching is exact 
in the LcrM in the MFA. 

Using the similar argument given in Sec lIVB] one can 
also show the phase quenching in the conventional LcrM 
with the SV{Nc) X SV{Nf)L x SU(A^/)fl x U(1)b sym- 
metries. The Lagrangian is given by 



where 



and the projection for mesons is to throw away d\ and 
from SO(27Vc) LctMb- It is easy to see that this pro- 
jection maps S0(2A^c) LcrMs into SU(A^c) LcrMg. 

Let us now consider how the orbifold equivalence can 
be shown within this model. The orbifold equivalence of 
the fermionic part is simple: because the fermionic part 
of the LcrM can be regarded as fermions in the presence 
of a background field cta, tta, d'^, and d^, the large- A^c 
equivalence holds as long as the background field does not ^^^^ 



(53) 



^B^^ [{d^^"'? + [d^of] + C/(a, ^'^), (54) 



(55) 



We note that our arguments below do not rely on the ansatz for 
A and /, as long as they are independent of Nc- e.g., A(T) = 
Ao [1 - (T/To)2] with some constants Ao = 0{N°) and Tq = 
0{N^) adopted in Ref. [48| . Our normalization is related to that 



in Ref. [481 Via. C ours ' 



U{a, n^) = ^ + i^"")' - iNcff] Ha, (56) 

where 7r° ~ ■ipi'y^T°'ip with r" being the SU(A^/) gener- 
ators. In the MFA, the effective potential is given by 
a summation of diagrams that have one central fermion 
loop with a number of meson fields a and tt" attached. 
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It is again only this fermion loop which depends on the 
chemical potential; the contributions of /i„ and /i^ are de- 
coupled, and the phase quenching is exact in the MFA. 



outside the pion condensation phase in the MFAI^ the 
phase quenching is exact for the free energy, chiral con- 
densate, and Polyakov loop in this model. 



D. Polyakov— Nambu—Jona-Lasinio model 



E. Polyakov-quark-meson model 



The Polyakov-Nambu-Jona-Lasinio (PNJL) model 
is an extended version of the NJL model by 
adding Polyakov-loop degrees of freedom to account for 
the confinement/deconfinement. The Polyakov loop (ex- 
pectation value) is defined by 



^(TrLt). 



Here L is an N^. x color matrix 



L(x) — Vexp 



i / (irA4(x, r) 
Jo 



(57) 



(58) 



with V being the path ordering, A4 — iAo, and /? = 1/T . 
The Lagrangian of the PNJL model is given by [l3 . [l6j 



^PNJL — -Ckin + -Cint + -^pot, 



(59) 
(60) 



where the interaction term Ci^t is taken to be the same as 
that of the NJL model, e.g., Eq. (pij) . On the other hand, 
84 in the kinetic term of the NJL model is replaced by the 
covariant derivative D4 = 84 — lA^ in (other deriva- 
tives are untouched, Di — di) and the Polyakov loop 
potential £pot = 'Cpot(^, 7") is introduced. The pa- 
rameters of the Polyakov-loop potential are determined 
by fitting lattice simulation data at /i = and finite T, 
but the detailed form of the potential is irrelevant in this 
paper. In order for the Polyakov loop to take the same 
expectation value between S0(2A^c) theory and QCD at 
any T (at /i = 0) as required by the large- A'^c equivalence 
in Sec. IIII Al we assume the same potential between the 
two theories. (This should be so at the level of MFA 
as well as one-meson-loop corrections, see the remark 
on 1/A^c corrections to gluonic operators at the end of 
Sec. IIII Al ) This is a necessary input for the model to be 
consistent with the underlying gauge theories. Once this 
assumption is made, the equivalence in the gauge sector 
of the model is trivially valid. 

On the other hand, the fermionic part of the PNJL 
model can be regarded as the NJL model in the presence 
of a background field A4 . As long as the background field 
does not break the projection symmetry (and it must 
be so because the Polyakov-loop potential is chosen at 
/I = where the projection symmetry is not broken), the 
equivalence in the fermionic sector is also satisfied from 
the argument in Sec. IIVBI 

Actually, as noted in Ref. [s^], the effective potential 
of the PNJL model satisfies the relation 



PNJL 



(M-b)U7=0 — i^PNJL(M/)UB=0, 



(61) 



One can also consider the extended version of the lin- 
ear sigma model by taking into account Polyakov loop 
degrees of freedom. This is known as the Polyakov-quark- 
meson (PQM) model [l3|. The Lagrangian is 



L 



PQM 



= c 



LcrM 



L 



pot J 



(62) 



where Ci^aM is the Lagrangian given in Eq. (|44|) with 
84 replaced by the covariant derivative D4 — 84 — iA4 
and Cpot is the same potential of Polyakov loop used 
in Eq. (|59l) . The proof of the LcrM can be extended 
straightforwardly, just as the proof of the NJL can be 
extended to that of PNJL. 



F. Chiral random matrix model 

In this subsection, we explain the orbifold equivalence 
of the chiral random matrix model (xRMM) following 
Ref. 01. The partition function of the xRMM 18] (for a 
review, see, e.g., Ref. [5^) is given by an integral over a 
Gaussian random matrix ensemble. 



Z = d^Y[detVf e 



(63) 



where $ is an x A^ random matrix element. The pa- 
rameter G is a normalization of the Gaussian. This the- 
ory does not have spacetime dependence; the size of the 
matrix A" corresponds to the spacetime volume which is 
taken infinity in the end (thermodynamic limit). 

The matrix structure of the Dirac operator T) is chosen 
such that it reproduces correct anti-unitary symmetries 
and global symmetry breaking pattern of the system. We 
can also add the quark mass m. quark chemical potential 
/i 0, and temperature T [Ullil into V. At T = 0, the 
Dirac operator is written as 



mfl ^ + ^fl 



(64) 



Here the matrix $ is taken to be real, complex or quater- 
nion real. Each case is respectively characterized by the 



There is an ambiguity to take the MFA of the PNJL model in 
the literature: its effective potential is complex at fj, 0, and 
thus, £ and £ are independent in general. However, the orbifold 
equivalence in the underlying QCD predicts that the Polyakov 
loop in QCDg agree with that in QCD/ outside the pion con- 
densation phase where 1 = 1. Therefore, the correct MFA in the 
PNJL model must satisfy 1 = 1 in QCDg in that region. This is 
consistent with the claim in Ref. 1511. 
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Dyson index /3 = l,/3 = 2or/3 = 4, which represents in- 
dependent degrees of freedom per each matrix element. 
^ = 1 corresponds to SU(2) QCD and Sp(2A^c) gauge 
theory, ^ = 2 to QCD with A^^ > 3, /3 = 4 to QCD with 
adjoint fermions and S0(2A'^c) gauge theory. 

Because the xRMM is a large- (not large- A^c) matrix 
model, one can prove the orbifold equivalence as in the 
same way as the field theories, just by replacing with 
N . In the following, we construct the orbifold projection 
from ^ = 4 RMM at finite (/? = 4 RMMs) to /? = 2 
RMM at finite = 2 RMM^) or ^ = 2 RMM at 

finite /i/ (/3 = 2 RMM/) at T = 0, which can easily be 
generalized to nonzero T (for the orbifold projection from 
/3 = 1 to ^ = 2, see Ref. [3].) The construction of the 
orbifold projections is almost the same as the projections 
from SO(27Vc)b to QCDb or QCD/. For simplicity, we 
consider degenerate quark masses mf — m. 

The partition function of the /3 = 4 RMM b is given by 



Z = / d^d^- e" 



S = Sf 



with 

and 

Sp = 



5„^^G^Tr$tci,, 



(65) 



(66) 



Nf 

E 

/=i 



V 



ml2N ^ + ^J''i-2N 



(67) 



where $ is a 2N x 2N quaternion real matrix and are 
complex 4A'^-component fermions. They can be decom- 
posed as 



effective potential of two-flavor [3 = 2 RMM with the 
quark chemical potentials and fid is computed, using 
the saddle point approximation for — > oo as [231 • 

f^RMM = G^[(o-« - muf + i<7d - mdY + 2{p - Xf] 

^ ln[(cr„ + fiu± iT){ad - T iT) + p^] 



: [{<Ju - fJ-uT iT){ad + fJ-d± iT) 



(72) 



The chiral condensate and pion condensate are related to 
au,d and p as 



2N " 



-G (Ju, 



-dx\nZ 



'G^p. 



(73) 
(74) 



A=0 



Outside the pion condensation phase, p = 0, the potential 
(l72t satisfies the relation: 



(75) 



By differentiating with respect to the quark mass, it fol- 
lows that the chiral condensates are identical between 
RMM/3 and RMM/ for p = 0. 



G. Strong-coupling lattice QCD 

The orbifold equivalence can be extended to the strong- 
coupling expansion of lattice QCD. For clarity, consider 
the action on the lattice with staggered fermions in the 
chiral limit [19]: 




aP + ia^ + ia} 
— + ia} aP — ia? 




(68) 



(69) 



Here "iAfl.L are 2A^-component fermions which are further 
decomposed into A'^-component fermions S,r^l and Cr,Li 
and a^ are N x N real matrices. 

In order to obtain /3 = 2 RMMs, we impose the pro- 
jection condition as 



(70) 



where J = —ia2 ® Iat and u> = e''^/^ as defined before. 
Then it is easy to see (3 = 2 RMM b is obtained after the 
projection. In the same way, (3 = 2 RMM/ is obtained 
by using 



(71) 



where Ji acts on the flavor indices. 

Let us check the exactness of the phase quenching in 
this model which is already observed in Ref. [28|. The 



S[U,x,x] = SG[U]+SF[U,x.xh 



1 ~ ^ReTr Uf,,{x) 



Here 



(76) 



Sf[U,x,x] 

= ^ E ^o(^) [x{x)e''Uo{x)x{x + 6) 

X 

1 1 

X j = l 

x{x)U,{x)xix+j) - xix + 3)uUx)xix)\ . (77) 



U,,,{x) = Ul{x)Ul{x + v)U,{x + p)U^{x), (78) 



is the plaquette, Ufj, is the SU(A^c) gauge hnk variable, 
X is the fermion field, 77p(x) is defined as rjo^x) = 1 and 

rij{x) = (— l)^i=i ^'-1 , and d is the number of spatial 
directions. 
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In the strong-coupling limit 5 — )■ oo, the gluon action 
Sg can be dropped, and the theory is given just by Sp. 
Because the orbifold projection can be defined for the re- 
sultant action at the lattice level, the equivalence imme- 
diately follows in the large- A^c limit. Note here that the 
large- iVc limit is taken for the action S after the strong- 
coupling limit g — ^ 00. In this case, however, the leading 
order in 1/Nc expansion does not correspond to the MFA 
in the literature j55| . 

The phase quenching is also exact at large to the 
next-to-leading order (NLO) and higher order in 
We again note that the large- A^c limit is taken after we 
truncate into the NLO (or higher) action of the strong- 
coupling lattice QCD. (For the higher-order calculation 
in l/g^, see, e.g., Ref. |56l|.) 

V. EQUIVALENCES IN HOLOGRAPHIC 
MODELS OF QCD 

In this section, we apply the orbifold equivalence to 
holographic models of QCD. Since the holography (or 
the gauge/gravity duality) maps a four-dimensional (4D) 
strongly-coupled gauge theory to a five dimensional clas- 
sical gravity theory in the large Nc and large 't Hooft 
coupling limits, we expect that we can use the large- iVc 
orbifold equivalence in holographic models, as originally 
proposed in Ref. Q- We consider below are the D3/D7 
model and the Sakai-Sugimoto model. 

A. D3/D7 model 

In this section we explain the equivalence in the D3/D7 
model ^ following Ref. [H. Let us start with 4D TV = 4 
\J{2Nc) supersymmetric Yang-Mills theory, which is re- 
alized around a stack of 2Nc D3-branes. The massless 
spectrum of D3-branes involves a vector multiplet on the 
worldvolume A0123 and three complex scalar multiplets 
describing the transverse motion X45, Xgy, X^g- At large 
Nc, this system is described by the type IIB superstring 
on AdSs X S^. In order to introduce 2Nf fundamental 
matters, we add 2Nf D7-branes, which wrap on C 
[Ifli: 
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Then two 2Nc x 2Nf chiral multiplets H describing 
strings from D3 to D7-branes and the reversed strings 
Ha = ^abH^'^ emerge. In the large- A'^j, limit with fixed 
Nf, one can neglect the back reaction and treat the D7- 
branes as probes in AdS^ x background^ By writing 



the AdSc, X metric as 

ds^ = ^S'^t.^dx^^dx'^ + E ^2^'' (79) 

the D7 are localized at yg = j/g = and extend along all 
the other directions. Then open strings connecting D3 
and D7 provide us with \]{Nc) x \}{Nf ) bi-fundamental 
matters, which resemble the U(7Vc) fundamental matters 
with U(7V/) flavor symmetry. The dynamics of quarks 
and mesons is described by the Dirac-Born-Infcld (DBI) 
action 

S-DBi = -Tt j d^S, TrV-det {G + 2iTa'F), (80) 

where ^ are the world-volume coordinates, G is the pull- 
back of the spacetime metric to the world volume and F 
is the field strength of the gauge fields on the brane, and 
T-j is the D7-brane tension. The chemical potential can 
be introduced as a background field of zeroth component 
of the gauge field on D7-branes. Here we choose the 
isospin chemical potential, 

* background ■ 7 /oi\ 

= =z;iJ2^^, (81) 

where J2Nf — —io'2'i-Nj- Starting with this theory, one 
can obtain an SO(27Vc) theory with Nf flavors at finite 
fiB via an orientifold projection, and a U(A^c) theory with 
Nf flavors at finite /is by further orbifold projection. For 
the orientifold projection, we introduce an 07-plane and 
a Z2 singularity as follows: 
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The geometric effect of the Z2 action is a reflection in 
the transverse directions a;6.7,8,9^ Hence the orientifold 
projection for the fields on D3 is 

^0123 ^ ^ (^0123)"'" ' 
^67,89 — (^67,89)"^ ■ 

Therefore, the orientifold projection for the gauge field is 

A'^^^{A,~Aj:), (83) 

so the projected gauge field is antisymmetric and spans 
an S0(2A^c) algebra. The field X45 is in an antisymmet- 
ric (adjoint) representation, while for the fields ^"67,89 the 



Note that the probe approximation in holography is different 
from the quenched approximation in QCD in general; while 



the latter does not distinguish /is and /ij as we have seen in 
Sec, nil Cl the former does distinguish them as we will see below. 
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orientifold action projects them to a symmetric represen- 
tation. Open strings connecting D3 and D7 are projected 
as 



and the fields on D7 are projected as 

-J2Nf (^0123)"^ '^27Vj' 



A' — 
^0123 — 



A' 



^T j-i 



^6789 — —J2Nf (^6789/' 

The chemical potential remains unchanged, 



J2N, 



(84) 

(85) 
(86) 
(87) 

(88) 



It can be regarded as both and /x/, because there is 
no difference between the two in the S0(2A^c) theory. 

By further performing a Z2 orbifold projection, one 
obtains a \]{Nc) theory with Nj flavors at finite ^b- The 
Z2 projection is 



A" 

^0123 
^67,89 



J2N^ ^0123 ^^2^' 
J2N^ X'l^ , 



— —J2Na ^1 



Jo 



67,89 "^2^^' 



for D3-D3 strings, 
for D3-D7 strings and 



jIlA T-l 



A" — T A" 
^0123 ~ •-'2Nj ^0123 ''2A'/' 



X'lz — J2Nf X'l^ J2Nf^ 



A" 

^6789 



-J2Nf ^6789 "^2^/' 



(89) 

(90) 

(91) 
(92) 
(93) 



for D7-D7 strings. The background field turns into the 
one corresponding to ^lb- The dual gravity geometry 
changes to AdS^ x RP^ through these projections. D7's, 
which were wrapping on C before the projections, 
wrap on MP^ ^ j^p5^ Because S'^ and RP^ are lo- 
cally the same, the DBI actions are the same except that 
the integration is performed on MP"^ instead of and 
the gauge fields are restricted to satisfy the projection 
conditionl^ Therefore the equations of motion derived 
from the DBI actions are the same unless the solutions in 
\]{2Nc) and S0(2A'c) theories break the projection sym- 
metry (or equivalently, if the solution is invariant under 
the projection symmetry), and hence the large- A^c equiv- 
alence holds. 
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FIG. 7. Phase diagram of the D3/D7 model (a) with (xb [53] 
and (b) with [5a - l6(]| | with a nonzero quark mass. The 
equivalence holds outside the rho meson condensation phase. 



The phase diagrams of the D3/D7 models with ^b |57| 
and /i/[58'-^6^ have been studied^l Schematic picture of 
the phase diagrams with a nonzero quark mass is shown 
in Fig. [T] There are phases with no charge density (anal- 
ogous to gluon plasma and a gas of mesons) and with 
nonzero charge density (quark-gluon plasma). With /i/, 
there is yet another phase, a rho meson condensation 
where U(l)/ symmetry is broken (analogous to the pion 
condensation phase in QCD/). In that region the equiv- 
alence does not hold. 



B. Sakai-Sugimoto model 

Sakai-Sugimoto model [2lj has reproduced the low- 
energy hadron spectrum successfully. It has also been 
used to study the chiral phase transition (see, e.g., 
Refs. j6ll - l63j ). It consists of Nc D4-branes wrapping 



As is well-known, the DBI action has an ambiguity of the order- 
ing of matrix variables. Here we assume there is a right ordering 
(though we do not know it explicitly) and the projections do not 
affect that ordering. 



Notice that the SU(Afc) theory with we are considering, which 
is dual to the string theory on AdS^ X RP^, is slightly different 
from the one studied in Ref . [57| , which contains three complex 
adjoint scalars and is dual to the string theory on AdS^ X S"". 
Ifowever these theories are equivalent in the large- A^c limit and 
the solution to the classical equations of motion are the same. 
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on a compactified circle, Nf D8-branes and Nf anti- 
D8-branes. Gauge symmetries on D8 and anti-D8 are 
identified as flavor symmetries U(A^/)l and \J{Nf)fi, re- 
spectively. When Nf/Nc ^ 1, D8 and anti-D8 can be 
treated as probes on the D4 background. In this setup 
it has been shown that a D8-brane and an anti-D8-brane 
merge to form single D8-brane, so that Nf D8-branes re- 
main and \J{Nf)L x \J{Nf)]i is broken down to U(7V/)v. 
This is the geometric realization of the spontaneous chi- 
ral symmetry breakdown. The flavor dynamics such as 
the meson spectrum can be read off from the DBI action 
for D8-branes and the Chern-Simons term. 

As in the D3/D7 model, we consider the probe (anti- 
)D8-branes in the D4-brane background. In the chiral 
symmetry breaking phase, the D4-brane background is 
given by 



\fiU) 



(94) 



where f{U) = 1 — U^^/U^. Here Ukk is a constant 
which is related to the radius of the compactified dimen- 
sion. D4-branes are extended along a;^ and a;^ directions, 
where x^ with /i = 0,1,2,3 is 4 dimensional spacetime 
and direction is compactified to . D8-branes are em- 
bedded in x^ , z, and S"' directions, where ^-direction is 
one dimensional space embedded in {U, x^) space. Extra- 
dimensions of are integrated out, and gauge fields on 
D8-branes and are related to (axial-)vector mesons 
and a pseudoscalar pion, respectively. Here, we explain 
the orbifold equivalence in this setup. It is straightfor- 
ward to extend the equivalence to other applications of 
the Sakai-Sugimoto model, for example, analysis of chi- 
ral phase structure, as long as it does not depend on the 
details of extra-dimensions. 

The equivalence can be shown in a similar fashion to 
the D3/D7 model. Normalizable modes of the gauge 
fields on the D8-branes correspond to the quark cur- 
rent and non-normalizable modes give their source. The 
chemical potential can be introduced as a non-vanishing 
background for the time component of the gauge fields. 
The background is taken to be proportional to Inj for 
^B, and proportional to a2 ® for [s^, [S^- We 

start with the Sakai-Sugimoto model with /i/. By per- 
forming the orientifold projection we obtain 0(2iVc) ana- 
logue of the Sakai-Sugimoto model with /i^. By further 
imposing the orbifold projection on the 0{2Nc) model, 
we obtain SU(iVc) model with /i^. Since /is is not com- 
patible with the symmetry for the orientifold projection, 
we cannot start with the Sakai-Sugimoto model with . 

The orbifold projection g, 



(95) 



acts on the D8 gauge fields as 

A,{x'^,x^) ^jig)A^ix^,-x^h-\g), 
A,(x^x^) = -j(g)A,ix^,-x^h~\g), (96) 

and similarly for the scalars. By using the gauge symme- 
try, 7((?) can be taken as 7(17) — (T2 (8) Itv/ ■ For the orien- 
tifold projection, the worldsheet reflection is taken in ad- 
dition to the spacetime reflection. It takes the transpose 
of the Chan-Paton factors and gives an additional sign 
for the gauge fields (but no additional sign for scalars). 
The orientifold projection acts as 



A^{x^ , 
Ai{x^ , 



^j{n)Al{x'^,-x^h-\n), 

7(^^)Af(x^-a;^)7-l(^^)• 



(97) 



There are two options for the orientifold, 7+(r2) = ^2Nf 
and 7-(f2) = J2Nf up to the gauge transformation, which 
give 0(2iVc) theory with 0(2iV/) flavor symmetry and 
Sp(2Afc) theory with Sp{2Nf) flavor symmetry, respec- 
tively. 

The orientifold projection of Sakai-Sugimoto model is 
studied in Ref. [6J]. In this case, we construct 0{2Nc) 
model. It can be obtained by introducing 06^ planes as 
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Then the orientifold projection acts on the gauge fields 
on the D8-branes as 



A^{x^,z) = -7+(^7)v4^(x^-z)7T^(^), 

^,(x^z) = 7+(^^)AJ(x^-z)7;l(^!). 



(98) 



Next we impose the orbifold projection. The fixed 
plane lies in the following directions: 
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The gauge fields on the D4-brane after the compactifi- 
cation become those for \J{Nc) symmetry. The orbifold 
acts on the gauge fields on the D8-brane as 

A^(x^ x') = 7(Z2)A^(a;^ -a;*)7-^(Z2), 
A,(x^ x') - 7(Z2)A,(x^ -x*)7-i(Z2), (99) 

where i = 5, 6, 7. 

We focus on the constant modes on S"'. By taking the 
orbifold projection first, U(27V/) gauge symmetry on the 
D8-brane is broken to U(A^/) x U(iVy), and we obtain two 
gauge fields A-^ and A"^ for each \J{Nf). The orientifold 
projection imposes the relation for these gauge fields as 
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Al{xf',z) = -A2(a;^,-z) and AKx^^^z) = Al{x^',-z). 
By imposing the both projections, we obtain the same 
effective theory with half flavors, but /i/ becomes /i^. 
Therefore, the equivalence holds as long as we consider 
only Z2 invariant sectors. 



VI. NUMERICAL EVIDENCE OF THE PHASE 
QUENCHING 

In this section we look at previous numerical simu- 
lations which compared QCD^ with QCD/. We shall 
confirm that QCD at large Nc and model calculations in 
the MFA provide us with a good approximation for the 
phase quenching in three-color QCD. 



A. Reweighting 

In Ref. 65], QCDs and QCD/ are studied by using 
the canonical formalism as a function of the number of 
up quarks, Q. The result of the former is obtained by 
the reweighting method [13, [IS] • They use two staggered 
fermions (corresponding to degenerate four up and four 
down quark species) with the bare quark mass am = 
0.14 on a 8"^ X 4 lattice. The canonical partition function 
Zc(Q) is obtained from the grand canonical partition 
function Zqc{iii) via the fugacity expansion. 



ZGc{V,T,fii)^Y.^c{V,T,Q)f 



(100) 



where 



ZGciV,T,fii)^ J dAdetD{fii)e-^^", (101) 
ZciV,T,Q)= [ dAdctQe-^™, (102) 



with detg being the projected determinant for the fixed 
quark number Q. From the above relations, the quan- 
tities detg and detD(/i/) are assumed to be related 
through 



det D{^j 



E 

Q 



detoe^'^^/^. 



(103) 



This relation allows us to extract detg, and hence, 
ZciQ)- The canonical free energy is then given by 
FcIq) = -{l/T)\nZc{Q). The canonical partition 
function and free energy at finite fiB can be obtained 
in a similar way. 

In the right panel of Fig. 1 and the left panel of Fig. 4 of 
Ref. [65l | , the free energy at various temperatures between 
0.5Tc and l.lTc are plotted as functions of Q. By putting 
these plots on top of each other, one can see a very nice 
agreement near the critical temperature and Q < 100. 
It clearly shows the validity of the phase quenching. It 
should also be remarked that the corrections are still tiny 



for Nf = 8, a larger number of flavors than Nf =2-1- 
1 in the real world [remember that the corrections are 
0{Nf/Nc) from our large- A^c argument in Sec. lIII Aj . 

In Ref. [66], three-color and two- flavor QCD^ and 
QCD/ are studied using staggered fermions with the bare 
quark mass am = 0.05 on a 8'^ x 4 lattice. The former is 
obtained by the phase reweighting from the latter. The 
chiral condensate and the Polyakov loop are computed 
for afi = 0.1 and a/z = 0.2, and the results of QCD^ and 
QCD/ agree within numerical errors, even for the average 
phase factor ~ 0.7. 



B. Imaginary chemical potential method 

The sign problem is absent when the chemical poten- 
tial is pure imaginary, fj, = iiiimg (Mimg G K) [l, [1^. This 
fact can be easily realized by an argument similar to the 
one around Eq. (jU; since the operator -|- i^img7^ 

is anti-Hermitian, its eigenvalues ±iA„ are pure imagi- 
nary, A„ € R, and the measure is positive semi-definite. 
Although the imaginary chemical potential is not phys- 
ical, it is useful if observables are analytic in fj? around 
/i^ = 0, because the values at /i^ > (real chemical 
potential), which are difficult to study due to the sign 
problem, may be obtained through an analytic continu- 
ation from < (imaginary chemical potential). Note 
however that the analyticity, which is necessary for the 
analytic continuation, can be lost at any phase transition, 
such as the chiral and deconfinement transitions. 

Our derivation for the large- A'^c equivalence in Sec. IIIII 
can also be applied for the imaginary baryon and 
isospin chemical potentials, {iJLu,Hd) — (*Mimg, *Mimg) and 
(Mu, Mrf) = (*Mimg, — «Mimg), without any modification. As 
a result, the chiral condensates (i/'V') b and {iptp)i take the 
same value at finite imaginary potentials as long as the 
projection symmetries are unbroken. 

In Ref. [67[, pseudo-critical temperatures of the chiral 
transition, T'c(/x), in two degenerate staggered fermions 
and three-color QCD at /x^ > were exploited by the 
extrapolations from /i^ < (for the bare mass am = 0.05 
on a 16'^ X 4 lattice). With an ansatz. 



Tc{0) 



l + ai 



they found j|i7,] 



ai 
ai 



-0.465(9) 
-0.515(11) 



for ^/, 
for 



(104) 



(105) 



which provide a nice quantitative agreement already at 
Nc = 3. As found from our arguments above, this differ- 
ence originates from the 1/Nc corrections. 



Roberge- Weiss periodicity 

At a finite imaginary baryon chemical potential, the 
jrand canonical partition function has the Roberge- Weiss 
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(RW) periodicity @ 

Z(^)=Z(^^ + ^) (neZ), (106) 

which can be understood as a generahzation of the center 
symmetry of the pure Yang-Mills theory; actually the 
Polyakov loop is transformed as £ — )■ e^'^™/^':^. 

In the confinement phase {£ = 0) the ground state 
also satisfies the RW periodicity. Therefore, in the large- 
Nc limit, there is no /ijmg-dependence, and thus, there 
is no /i-dependence at /x > until the phase transition 
happens. This is consistent with an important property 
of the large- QCD that observables are T-independent 
in the confinement phase[3 This can also be understood 
physically from the fact that 0{N^) observables, such 
as the chiral condensate, cannot be affected by thermal 
excitations of noninteracting mesons and glueballs, which 
have only 0(iV°) degrees of freedom (70| . 

In the deconfinement phase {£ ^ 0) the vacuum 
does not respect the RW periodicity and nontrivial /i- 
dependence can appear. 

C. Taylor expansion method 

Another common approach to circumvent the sign 
problem is the Taylor expansion method; one expands 
the expectation value of an observable in power series of 

fi/T [Big, 

(<^)B^±c^{^y (107) 

71=0 

in QCDs and 

oo 

in QCD/. Taylor coefficients and c^, which are func- 
tions of the temperature T, can be determined by the 
simulation at /x = 0. The large- iV^ equivalence tells that 
the coefficients agree in the large- A^c limit. 

In Ref. [26^1, the coefficients and for the chiral 
condensate and the pressure of the quark-gluon plasma 
have been calculated^ in three-color and two-flavor 



This can be understood as a generalization of the Eguchi-Kawai 
reduction [69| | stating that observables are independent of the 
size of the compactified direction in the confinement phase. If 
we use it for the compactified T direction, T-independence of 
observables immediately follows. It leads to the /^-independence 
at finite T for fi < fis/^c, because there is no /i-dependence at 
T = for /t < fis/Nc. 

For odd n, and vanish, and the first nontrivial /t- 
dependences appear in and Cj. Although (n > 4) have 
been calculated, (n > 4) have not been calculated in Ref. l26l . 
(Note that, for n > 4, they use the same symbol c^j for another 
quantity.) 



QCD. Their calculations are performed using staggered 
fermions with the bare quark mass am = 0.1 on a 16^ x 4 
lattice. The coefficients for the pressure are [2^ 



T/T, 


^2 




0.81 


0.0450(20) 


0.0874(8) 


0.90 


0.1015(24) 


0.1551(14) 


1.00 


0.3501(32) 


0.3822(26) 


1.07 


0.5824(23) 


0.5972(21) 


1.16 


0.7091(15) 


0.7156(14) 


1.36 


0.7880(11) 


0.7906(9) 


1.65 


0.8157(8) 


0.8169(7) 


1.98 


0.8230(7) 


0.8250(6) 



Although the difference between cf and C2 are not very 
small for T < Tc in the chiral symmetry broken (and 
confined) phase, they agree exceptionally well for T > 
Tc- This tendency can naturally be understood, as we 
have argued in the end of Sec. IIII Al The coefficients for 
the chiral condensate are shown in the second panel of 
Fig. 3.6 of Ref. [26]. There the agreement is even better; 
the coefficients agree within errors for T/Tc > 0.87. 

VII. DISCUSSIONS AND OUTLOOK 

In this paper, we have systematically developed the 
string-inspired technique of the orbifold equivalence in 
effective models, holographic models, and lattice meth- 
ods for QCD. As a consequence, we provided the criteria 
for the validity of the phase quenching, as summarized in 
Sec. m The phase quenching does not produce any quan- 
titative difference of the chiral and deconfinement phase 
transitions in the mean-field approximation (MFA) and 
to the one-meson-loop corrections, respectively, outside 
the pion condensation phase of the phase-quenched the- 
ory. 

In the pion condensation phase, the orbifold equiva- 
lence breaks down. Also the 1/Nc expansion itself may 
no longer capture the physics of real QCD[^ Since the 
validity of the MFA in model calculations may be tightly 
connected to the validity of the 1/Nc expansion as we 
have seen in Sec. IIV Al it is possible that the 1/A^c ex- 
pansion as well as the MFA are not useful inside the pion 
condensation phase of QCD/. If so, previous model cal- 
culations supporting the existence of the QCD critical 
point (see Ref. [76] for a review) may not be reliable, be- 
cause it was observed only inside the pion condensation 



In the large-TVc limit, there is no nuclear liquid-gas transition 
[71II and no color superconductivity [72| which are expected to 
be realized in real QCD. This necessitates a phase transition as 
a function of Nc from homogeneous to inhomogencous matter 
(from a nuclear gas or nuclear liquid to a nuclear crystal 73, 7^ 
and from a color superconductor to a chiral density wave 72j or 
a chiral quarkyonic spiral |75^|1. 
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phase in model calculations under the MFA [36|. Actu- 
ally, by utilizing the large- iVc equivalence, it has recently 
been shown that QCD critical point cannot exist outside 
the pion condensation phase in the large-iVc QCD and 
effective models in the MFA [tJ- Therefore, the effects 
beyond the leading order in l/Nc or those beyond the 
MFA should be taken into account to describe the real- 
istic dense matter. 

Of course it is important to understand the fate of 
the chiral phase transition outside the pion condensation 
phase. For example, one could study the curvature of the 
chiral critical surface [ts] away from ^ ^ 0, from which 
we may hopefully infer the behavior of the chiral phase 
transition and the (non) existence of the possible QCD 
critical point at larger /i. Based on the phase quenching 
approximation and using the rooted staggered fermions, 
it was numerically suggested in Ref. [73 that the QCD 
critical point does not exist in three-flavor QCD^ for 
/i ^ m.^/2. A more decisive conclusion should be drawn 
by detailed numerical calculations in the future. 

Note added. — After this work was completed, we were 
informed about Ref. (80| which claims to prove the non- 
perturbative equivalence between QCD^ and QCD/ at 
large Nc- However, their argument and their conclusion 
about the l/iVj, corrections are nothing more than the 
previous perturbative arguments (sol. l3ll|. They consider 
the expansion of the fermion determinant in 1 /Nc around 
the Yang-Mills vacuum without quarks (i.e., quenched 
approximation of QCD^), which is nonperturbatively 
different from QCD^ itself, as argued in Sec. IIII CI of 
our present paper. Note that nonperturbative fermionic 
condensates do appear in the large- iVc limit, as opposed 
to their claim, which falsifies their assumption. Note also 
that the existence of condensates itself does not spoil the 



equivalence. Rather the symmetry realization is the key; 
for example the equivalence holds with a nonzero chiral 
condensate. They also missed quantities for which the 
phase quenching becomes exact at large Nc- 

In order to give a nonperturbative proof of the equiva- 
lence, one has to take into account the symmetry realiza- 
tion, as emphasized in Ref. [Tlj (see also Refs. [sl. R [33| ) . 
Probably the mostproniising approach is to use the orb- 
ifold equivalence [3, |j| and to consider the Hamiltonian 
dynamics of the theory truncated at one-fermion-loop 
level, by closely following Ref. [ll|. We defer it to a 
future problem. 
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